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Abstract 


Starting from a simple empirical parametrization of the scattering amplitude, successfully de¬ 
scribing the dip-bump structure of elastic pp scattering in t at fixed values of s, we construct a toy 
model interpolating between missing energy intervals to extract the Odderon contribution from 
the difference between pp and pp elastic and total cross sections. The model is fitted to data from 
y/s = 23.5 GeV to 7 TeV and used to extract the Odderon and its ratio to the Pomeron. From our 
fits, a unit intercept Odderon follows, as predicted by J. Bartels, L.N. Lipatov, and G.P. Vacca, on 
the basis of perturbative quantum chromodynamics. 
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I. INTRODUCTION 


The nature of the Odderon — an asymptotic odd-U Regge pole exchange, counterpart 
of the Pomeron — for a long time remains a subject of debate. Although there is little 
doubt about its existence, we still lack direct evidence of the Odderon. Various reactions 
supposedly dominated by Odderon exchanges, called “Odderon hlters” [l|, may offer only 
indirect evidence either because of low statistics or contamination by competing exchanges. 

In quantum chromodynamics the Odderon corresponds to the exchange of an odd num¬ 
ber of gluons. Relevant calculations were done in a number of papers; see and earlier 
references therein. 

The only direct way to see the Odderon is by comparing particle and antiparticle scat¬ 
tering at high enough energies. The high-energy proton-proton and proton-antiproton elas¬ 
tic scattering amplitude is a difference or sum of even and odd U-parity contributions, 
t) = ’’'Even" ± '^Odd”, where, essentially, the even part consists of the Pomeron and / 
Reggeon, while the odd part contains the Odderon and the u Reggeon. It is clear from the 
above formula that the odd component of the amplitude can be extracted from the difference 
of the proton-antiproton and proton-proton scattering amplitudes, and, since at high enough 
energies the contributions from secondary Regge trajectories die out, this difference offers 
a direct way of extracting the Odderon contribution. Unfortunatelly, pp and pp elasctic 
scatterings were typically measured at different with the exception of the ISR energies 
of 31, 53, 62 GeV (see Fig. [I]). 

At present, the only way to extract the Odderon from the difference of the pp and pp 
scattering amplitudes is by means of a reliable interpolation of both amplitudes (or cross 
sections) over the missing energy regions. While the energy dependence of the forward 
amplitude (or total cross sections) is controlled by the Regge pole theory, its t dependence, 
especially in the dip-bump region, to a large extent is model dependent and unpredictable. 

A simple, general, and reliable parametrization of the complicated diffraction structure at 
high energies at any hxed energy is a sum of two exponentials in t related by a complex phase 
6*“^. Using this generic expression, Phillips and Barger (PB) (for brevity we shall refer 
to it as the PB ansatz) obtained good fits to the proton-proton differential cross sections, 
including the dip-bump region at several CERN ISR fixed energies. In Refs. 5l-l7|] the model 
was extended and improved, in particular, by accounting for the nonexponential behavior 
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da/dt data for pp, pp or both 



year 


FIG. 1. Timeline of proton and antiproton elastic scattering measurements. New accelerators are 
run first at the maximum available energies; however, at the start of the SppS accelerator, the pp 
and the pp elastic scattering data were measured at the same y/s = 31, 53 and 62 GeV. 


in the low |t| region. In Ref. 9| the parametrization was shown to be applicable to proton- 
antiproton scattering as well. Its connection with inelastic reactions is discussed on p. 185 


of Ref. 


10 |. 


The PB ansatz reads 


^(s, t) = i[\/]4exp(Rf/2) exp(i0(s))\/Gexp(T)t/2)], 


( 1 ) 


where s and t are the standard Mandelstam variables; and 0 were fitted to each 

energy independently; i.e., energy dependence in the PB ansatz enters parametrically. 

The success of the simple PB parametrization motivates its further improvement, exten¬ 
sion, and utilization. In papers js, 6| the low-|f| behavior of the PB ansatz was improved by 
modifying its simple exponential behavior by (a) inclusion of a two-pion threshold required 
by analyticity and (b) by means of a multmlicative factor reflecting the proton form factor. 
Achieving good hts to the TOTEM data j^, at the LHC energies of y/s = 7 TeV, the mod- 
ihed ansatz was used to predict the behavior of the observables at future energies as well as 
the expected asymptotic behavior of the cross sections. 
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The model was also tested in j?, 9| against the TOTEM data 4], and, contrary to many 
alternative models, it works reasonably well. In Ref. jh] the PB ansatz Eq. o was improved; 
in particular, its small-t behavior was modified. The fitted parameters are collected in 
Table II of that paper. For example, the authors of jb] quote the following values for these 


parameters for pp scattering at ^/s = 53 GeV: vCd = 6.55, \/C = 0.034 in y mb/GeV^ and 
B = 10.20, D = 1.7 in GeV“^, 0 = 2.53 rad. 

As already mentioned, in addition to elastic pp scattering, the PB ansatz also describes 
pp data, with a different set of the parameters (see below), thus opening the way to be used 
as a tool in extracting the Odderon from the difference of the two. However, in its original 
form, the PB ansatz does not describe the y/s dependence of the model parameters. 

In this paper, we try to remedy this limitation by combining the appealingly simple and 
efficient form of its t dependence with energy dependence inspired by the Regge pole model. 
Work in this direction was started in papers 

We address the following issues: (1) we smoothly interpolate between the values of the 
parameters fitted at fixed energy values, (2) extract the Odderon contribution from the 
difference of the pp and pp cross sections, and (3) compare the energy dependence of this 
difference with the prediction 12| based on perturbative quantum chromodynamics, by 
which the intercept of the Odderon trajectory equal to one. 

The structure of this paper is as follows: in the next section, we introduce the s/s 
dependence of the parameters of the PB model, which in the next section are determined 
from the fits to data. Then we discuss the results, in particular the s/s dependence of the 
results, including, for example, the Odderon contribution and the Odderon/Pomeron ratio. 
Finally, we summarize and conclude. 


II. THE GENERALIZED PB MODEL 

We use the norm where 

atot = 47rAA(f = 0) = 47r[\/A + \/C cos0] (2) 

and 

^ = 7r\A{t)\‘^ = 7r[Ae^^ + cos0]. (3) 

CLL 
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Following the Regge pole theory, we make the following assignment, 


\/]4 —)■ \/ R(s) = ais + 025 ^“^^ ^/c \/C{s) = (4) 

inspired by the Donnachie and Landshoff model of cross sections [see Eq. ([2])] with 
effective falling (subleading Reggeons) and rising (Pomeron) components. It follows from 
our hts that the falling (subleading Reggeon) components in \/C are small; hence, they are 
neglected. 

The slopes B and D in the Regge pole theory are unambiguously logarithmic in s, pro¬ 
viding shrinkage of the cone: 

B ^ B{s) = bo + hi ln(s/so), D D{s) = do + di ln(s/so). (5) 

In the above formulae a normalization factor sq = 1 GeV^ is implied. 

The phase (j) is the weakest point of this “toy” model or generalized PB model. In Regge 
theory, it should depend on t rather than on s. Fortunately, at high i/i the dependence of 0 
on energy is weak (see the hts below). However, this is not the case as the energy decreases. 
The best we can do is to ht the data with 


cos(0(s)) = ko + kis 


( 6 ) 


The “low”-energy behavior is a weak point in any case. Apart from the varying phase, we 
must account in some way for the subleading (/ and u) Reggeon contributions. This is 
done partly by the inclusion in of a decreasing term (absent in \fC). A complete 

treatment of these terms with proper t-dependent signatures will require a radical revision 
of the model, and we hope to come back to this issue in the future. 

Now we proceed with this simple approach that has a chance to be viable at high ener- 


13| and the above complications may be 


gies, where the Pomeron and Odderon dominate 
insignihcant. 

To understand better the existence of any connection between the ansatz ([1]) and the 
Regge pole model, we plot the values of the parameters A, R, G, D, and 0 against s and ht 
their “experimental” values to Regge-like formulas. 

This can be done in two complementary ways: A successive “two-step” ht. First, we ac¬ 
quire the values of the parameters A, R, G, R, 0 from the hts to the pp and pp data, then 
we ht their Regge forms (see below) to the obtained “experimental” values of A, R, G, R, (j). 
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Alternatively, one may determine the parameters of Eqs. (jl]) and ([5]) from a single simul¬ 
taneous fit to all available data. We chose the hrst option (5 parameters) since otherwise 
there were too many (at least 12) free parameters. Thus, we proceed with a two-step £t, by 
which the hnal values are determined from a £t to the “experimental” values of A, B, C, D, 
and <p. 

We fitted separately pp and pp in two variables, s and t, by using pp and pp data on total 
and differential cross sections ranging from the ISR to the LHC for pp and from SppS to the 
Tevatron for pp. 

Here the following remarks are in order: 

(1) It is clear from Eqs. (jl]) and ([5]) that, while the parameters A and C are particu¬ 
larly sensitive to the data on total cross sections, B and D are correlated mainly with the 
differential cross sections (the slopes). 

(2) Although we are interested mainly in the high-energy behavior (the Odderon), low- 
energy effects cannot be fully neglected. They are taken into account approximately by 
including in A and C subleading terms of the type e ~ 0.5. 

(3) Having htted A, B, C, D, and (j), we perform a cross-check by calculating the resulting 
total cross sections. 

(4) At high energies, the proton-proton and antiproton-proton total cross sections are 
supposed to converge. (We consider only this simple option, although we are aware of 
alternatives.) Since the existing data are not yet in this asymptotic domain, we introduce 
an extra constraint. 

(5) The most delicate issue is the phase, which in Re^e phenomenology is expected 
to be t, rather than s, dependent. Our hts (and those of |5|-l7|) show considerable energy 
dependence of the phase at low energies but weak dependence at high energies, where we are 
particularly interested in looking for the Odderon signal. Postponing the introduction of a 
true Regge-pole-motivated, t-dependent phase to a further study, here we assume a simple 
parametrization cos^ = fco + 


III. FITTING THE MODEL TO THE DATA 

We have calculated the s dependence of the parameters using a fitting strategy consisting 
of three consecutive steps described in detail below. In doing so, the following criteria were 
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applied: 


• best for oach fit; 

• the —t range was set within 0.35-2.5 GeV^; 

• for each hxed energy, the model was htted simultaneously to da/dt and a tot (focusing 
on the dip region); 

• from the resulting hts, atot was reconstructed in the whole available energy range. 


The values of the httec 
those obtained in Refs. |3|, 


arameters, in general, are consistent, within about 30%, with 
^ , except for the small value of -sfC deviating considerably from 


those quoted in the above papers, which used a different t range for their data analysis. 
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Energy 

(GeV) 

Va 

B 

Vc 

D 

cos{4>) 

xVndf 

23.4 

3.13± 0.6% 

8.66± 0.4% 

0.019± 8.3% 

1.54± 5.1% 

-0.97± 0.3% 

1.6 

30.5 

3.21± 0.2% 

8.95± 0.3% 

0.014± 7.4% 

1.28± 5.6% 

-0.98± 0.2% 

1.1 

44.6 

3.33± 0.7% 

9.32± 0.5% 

0.017± 8.0% 

1.45± 5.3% 

-0.93± 0.8% 

1.7 

52.8 

3.38± 0.3% 

9.44± 0.6% 

0.017± 7.6% 

1.43± 5.0% 

-0.92± 0.9% 

1.1 

62.0 

3.49± 0.5% 

9.66± 0.6% 

0.018± 9.9% 

1.53± 6.3% 

-0.92± 1.6% 

1.5 

7000.0 

8.51± 1.6% 

15.05± 0.8% 

0.670± 2.3% 

4.71± 0.8% 

-0.93± 0.3% 

1.4 


TABLE I. Values of the parameters from a fit to the pp data at various ^/s. The quoted errors 
correspond to the relative errors, as given by CERN MINUIT fitting package (status= converged, 
error matrix accurate). 


A. Step 1: Fitting the parameters to da/dt and atot for fixed energies 


Figure [2] shows a fit to the data on the pp and pp differential and total cross sections. The 
parameters A, B,C, D, and (j) were fitted to each energy separately. Given the simplicity of 
the model, the fits look reasonable. 


Figure [3] shows the fitted values of the parameters A, B, G, D, and (j) both for pp and pp 
scattering to be used as “experimental” data in the second stage of our fitting procedure, 
in which the explicit expressions (jl]), ([5]) and ([6]) are inserted. The fitted values of the 
parameters and relevant x^/^DF values are quoted in Tables HI and UTl 


Note the difference between the present model and that of Refs. 


BJzi 


in the low-|t| be¬ 


havior of the differential cross section, corrected in IS-Sl and taking account of the deviation 

n 

from an exponential. Normalization in Ref. [3[ is arbitrary since, in that paper, only the 
differential cross section is shown. Total cross sections are not shown in the above papers. 
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FIG. 2. The PB model fitted to the pp and pp data at discrete energy values. The da/dt data are 
from Refs. 


16l42ll| . The data on atot are from the Particle Data Group database 


22 |. 


Energy 

(GeV) 

Va 

B 

Vc 

D 

cos{4>) 

xVndf 

63 

3.43± 1.1% 

10.07± 1.3% 

0.022±30.8% 

1.90±14.8% 

-0.60±22.7% 

0.7 

546 

5.06± 1.2% 

11.25± 1.3% 

0.204±21.0% 

3.55± 8.6% 

-0.86± 2.7% 

0.6 

630 

5.13± 3.9% 

11.26± 3.7% 

0.176±26.6% 

3.23± 9.6% 

-0.81± 7.9% 

0.5 

1960 

6.85± 3.7% 

12.46± 3.3% 

0.629±41.6% 

4.69±15.4% 

-0.90± 3.6% 

0.4 


TABLE II. Values of the parameters fitted to pp data 


9 
























B. Step 2: Fitting the parameters a*, bi,... entering Eqs. Q, ([5]) and ([6]) to the 
“data” A, B, C, D, cos{4>)..., quoted in Tables I and II. 


The resulting values of the parameters after the second stage of fitting 


yjApp{s) = 

^JCpp{s) = 0.00117s°-3^®, 

Bpp{s) = 5.13 +0.5551ns, 
Dpp{s) = -0.838+ 0.312 Ins. 
cos(0pp(s)) = -0.928 - 0.863 s"°-^2^. 

^jApp{s) = 1.31s°'^°® + 4.28 s-°'2®® 
^JCpp{s) = 0.00177s°-^^®, 

Bpp{s) = 7.87 + 0.2741ns, 
Dpp{s) = -0.552+ 0.312 Ins, 
cos(0pp(s)) = -0.928 + 4.37 s“°-^2®. 


are 

( 7 ) 
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The lowest right icon in Fig. |3] is a “cross-check”, showing the pp and pp total cross 
sections calculated from Eq. ([2]) with the explicit values of the parameters defined by (jl|), 
()S|) and ([S]). The pp total cross section turns down at highest energies, deflecting dramatically 
from that of pp. The reason for that nonphysical effect is the scarcity of pp data, leaving 
too much freedom in the high-energy extrapolation of the cross section, where one expects 
asymptotic equality = crf^ for s ^ oo; see Eq. (47) in Ref. [1^. This deficiency should, 
and can, be cured by imposing an additional constraint on the model. This will be done in 
the next subsection, by fixing (tuning) the parameter (the leading powers in s of y/A) 
to be the same in pp and pp scattering. 
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FIG. 3. Energy-dependent values of the parameters extracted from a fit to pp and pp data. The 
data on atot are from the Particle Data Group database 22]. 
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C. Step 3: Tuning (refitting) the parameters by imposing the asymptotic con¬ 
straint 


't ) 


oo 


The above, unbiased fit does not satisfy automatically the required (see [1^) asymptotic 
constraint = cxf^ since the available freedom (especially due to the lack of simultaneous 
pp and pp elastic scattering data at \/s = 540, 630, 1800 and 7000 GeV) leaves much freedom 
for the extrapolation to energies beyond the existing accelerators. To remedy this problem, 
we have tuned the parameters to meet the above constraint in the currently available energy 
range. Below are the results of the “tuned” fit (see Fig. 0]) satisfying the asymptotic 
condition in the ^/s < 14 TeV energy range. 

The refitted s-dependent values of the parameters for pp and pp scatterings are 


^JApp{s) = 1 . 41 s °-° 9 ®® + 2 . 78 s -°- 2 ®^ ( 8 ) 

\JCpp{s) = 0.00223s°-^°®, 

Bpp(s) = 4.86 + 0.586 In s, 

Dpp{s) = -0.189 +0.2501ns. 
cos(0pp(s)) = -0.928 - 0.838s"°-^^^ 

^App{s) = 1.41s°-°^®® + 4.00 s"°-2®^ 

^JCpp{s) = 0.00588s°•2®^ 

Bpp{s) = 6.55+ 0.398 Ins, 

Dpp{s) = 2.351 +0.0681ns, 
cos{(f)pp{s)) = -0.908 + 4.376s"°-^2®. 
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FIG. 4. Energy-dependent values of the parameters from a fit to pp and pp data, constrained 
by as ^/s —>■ oo. The lowest right icon shows a cross-check for the (asymptotically 

converting) total cross sections. 
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Pomeron Odderon 





FIG. 5. Odd (Odderon) and even (Pomeron) parts of the da/dt cross sections calculated from the 
difference \A\^p — \A\pp = Aodd and from the sum \A\pp + \A\pp = T,pom fitted to the data. 

IV. THE ODDERON 


The existence of a parametrization for both pp and pp scattering offers the possibility to 
extract the odd-C* contribution by using the formula 


A^^ = Aeven ± A, 


PP 


odd") 


( 9 ) 


where Aeven and Aodd are, respectively, the G-even and G-odd components of the scattering 


amplitude; see Table 1 in Ref. 
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While the G-even component contains the Pomeron and the / trajectory (both known), 
the G-odd part is made of the poorly known Odderon and the familiar u trajectory. At 
the LHC energies, the contribution from secondary trajectories, e.g., / and u, is negligible 
(see Ref. jl^); therefore, by taking the difference between the known (fitted) pp and pp 
amplitudes, one gets a pure odd-G contribution that, in the LHC energy range, is the Odd- 
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Odderon 



FIG. 6. Real and imaginary parts of the difference calculated 

model fitted to the data. As seen from the figure, this difference tends to a 

□ 

corresponding to a unit intercept Odderon [1^ . 


from the present 
(small) constant, 


eron. From the explicit expressions for pp and pp amplitudes (cross sections), we calculate 
the Odderon amplitude (or its contribution to the cross section) by taking the difference 
j^pp _ j\pp = Aodd- The result (the energy dependence for several hxed values of t and t 
dependence for several hxed values of s) is shown in Fig. 4. 

The extracted model parameters are used then to evaluate the even and odd contributions 
to the forward scattering amplitude. In Fig. [5] we show the Pomeron and the Odderon 
contributions as the sum or the difference of the differential cross section of pp and pp elastic 
scattering. We see that, as expected, the Pomeron dominates at large colliding energies, 
while the Odderon contribution is small and at t = 0 even changes sign. A particularly 
interesting feature is shown on the lower right panel of Fig. |5l where the Odderon/Pomeron 
ratio is shown at different values of t at various y/s. Apparently, at y/s ~ 100 GeV, the 
Odderon/Pomeron ratio becomes t independent and the t-dependent curves pass through 
the same point of about 0.03. 

An important hnding of our paper is the near energy independence of the odd-G con¬ 
tribution to the scattering amplitude, which, at high (e.g., those of the LHC) energies, is 
dominated by the Odderon exchange. In Fig. [6] we plot the odd-G contribution calculated 
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from the difference ([9]). Both the real and imaginary parts of this difference at t = 0 tend 
to a (small) constant value, which correspond to a unit-intercept Odderon as predicted [12 1 
from quantum chromodynamics. 


V. CONCLUSIONS 


The present work is a semiquantitative estimate of the possible odd-C contribution to the 
scattering amplitude (cross sections) at high energies. “Semiquantitative” implies limitations 
due to the following: 


For the sake of simplicity, we ignore the low-t nonexponential behavior (sharpening) 
of the differential cross section. This simplihcation has dramatic impact on the low-t 
behavior of the Odderon contribution because of the large errors due to the cancellation 
of the Pomeron contribution. The importance of the low-t effects was emphasized, e.g., 
in Refs, js, 6|. 

The s- (rather than t-) dependent signature factor (phase) is in agreement with the 
data, but it is in contrast to expectations based on Regge phenomenology. 

The odd-U contribution to the amplitude, equal to the difference (|9]), shown in Fig. 
[6l tends to a constant limit determined by a unit-intercept Odderon, predicted by 
Bartels, Lipatov, and Vacca jl2 |. 


It is an oversimplihed treatment of the low-energy (“secondary Reggeons”) contribu¬ 
tions (Odderon and Pomeron here implies, generally, odd and even exchanges). 

There is the absence, for the moment, of any physical interpretation in terms of 
Reggeon exchanges of the components in the PB ansatz. 


Given these limitations and simplihcations, our approach can be considered as semiquan¬ 
titative, showing however some new aspects of the enigmatic Odderon. 

We strongly recommend to run the LHC accelerator at the injection energy y/s = 900 GeV 
and at the Tevatron energy of 1.8-1.96 TeV, so that the missing energy range of pp elastic 
scattering will be covered and elastic pp scattering data will be measured in the region where 
elastic pp scattering is already measured. It would also be desireable to measure elastic pp 
scattering at the \/s = 500 GeV region — the upper energy range of the RHIC accelerator. 
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Using the presently available data, the indehnite rise of the C(s), multiplied by a negative 
“signature factor” cos0, prevents the use of the generalized PB model beyond the LHC 
energy region. 

In the future we intend to 


rewrite the PB ansatz with correct, t-dependent Regge signature factors [h], remem¬ 
bering that in the present study the even and odd parts of the amplitude differ only 
by the values of the htted parameters, without the identihcation of particular Regge 
exchanges (trajectories); 


£t to the data to provide hints concerning the physical meaning of the components 
(two or more); for example, our hts (and those of Refs. 0, 5 -7]) indicate that A » C, 
but C rises with energy faster than A. Similar considerations may help in relating the 
second term with the Odderon or a “hard” Pomeron, for which the slopes {B and D) 
may be indicative. 


• add lower-lying contributions (subleading Reggeons) which are inevitable at lower 
energies, at the prise of giving up of the attractive simplicity of the model; 


use other asympto, 
the work in Ref. 


ic Pomeranchuk-like constraints in determining the parameters (see 
that points to this interesting direction); 


calculate and £t the model to other observables, e.g., the ratio of the real to imaginary 
parts of the amplitude p{s, t), of the slope B{s, t), etc. Its relation to inelastic processes 
(see 10|) may offer additional information. 
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